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O Abstract 

K/i We analyse the two-dimensional MHD configurations characterising the 

1 steady state of the accretion disk on a highly magnetised neutron star. The 
model we describe has a local character and represents the extension of the 

^ crystalline structure outlined in [1], dealing with a local model too, when a 

^ specific accretion rate is taken into account. We limit our attention to the lin- 

earised MHD formulation of the electromagnetic back-reaction characterising 
the equilibrium, by fixing the structure of the radial, vertical and azimuthal 
profiles. Since we deal with toroidal currents only, the consistency of the 
model is ensured by the presence of a small collisional effect, phenomenolog- 
ically described by a non-zero constant Nernst coefficient (thermal power of 
the plasma). Such an effect provides a proper balance of the electron force 
equation via non zero temperature gradients, related directly to the radial 
^ and vertical velocity components. 

^ We show that the obtained profile has the typical oscillating feature of 

the crystalline structure, reconciled with the presence of viscosity, associated 
to the differential rotation of the disk, and with a net accretion rate. In fact, 
we provide a direct relation between the electromagnetic reaction of the disk 
and the (no longer zero) increasing of its mass per unit time. The radial 
accretion component of the velocity results to be few orders of magnitude 
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below the equatorial sound velocity. Its oscillating-like character does not 
allow a real matter in-fall to the central object (an effect to be searched into 
non-linear MHD corrections), but it accounts for the out-coming of steady 
fluxes, favourable to the ring-like morphology of the disk. 

keywords: Accretion disks; Plasma physics; Nernst coefficient 

PACS Nos.: 97.10.Gz, 52.30.-q 

1 Introduction 

A long standing problem in astrophysics concerns the mechanism of accretion 
that compact objects manifest in the presence of lower dense companions. It is 
well established by observations [2] that such an accretion profile takes place often 
with the morphology of a thin disk configuration. Since from the very beginning of 
these studies, increasing interest raised about the details of the angular momentum 
transport across the disk structure. The solution of this problem appears settled 
down in the case of a compact object, for which the intrinsic magnetic field does not 
exert a significant Lorenz force on the elementary charged constituents of the in- 
falling plasma. In this hmit, the fluid-dynamics description appears well-grounded 
and the resulting axisymmetric configurations well understood (see among the first 
analyses of the problem[3, 4, 5]; for a relativistic analysis of plasma accreting into 
a black hole, see [6]). In the fluid-dynamics paradigm, the accretion mechanism 
relics on an angular momentum transfer which is allowed by the shear viscosity 
properties of the disk material. The differential angular rotation of different disk 
layers is associated with a non-zero viscosity coefficient, which accounts for the 
diffusion and turbulence phenomena emerging in the micro-scale structure of the 
fiuid. 

However, when the magnetic field of the central object is sufficiently high, the 
electromagnetic back-reaction of the disk plasma becomes very important. As 
shown in [1], the Lorenz force introduces a crucial couphng between the radial 
and the vertical equifibrium, which deeply alters the morphology of the system. 
In particular, the radial dependence of all the fundamental quantities acquires an 
oscillating character, modulating the background profile. Furthermore, as outlined 
in [7], when the disk plasma is characterised by a /3-parameter close to the unity, 
such an oscillating-like behaviour stands as the dominant effect and the disk is in- 
deed decomposed in a ring profile. However, such a two-MHD approach is pursued 
neglecting the accretion rate of the disk in the leading order and thus avoiding the 
problem concerning the azimuthal balance of the force acting on the electrons. In 
fact, in these two works, toroidal currents and matter fluxes arc addressed only. 

This question about a consistent electron force balance along the azimuthal 
direction (the azimuthal electric held is request to vanish by the axial symmetry) 
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affects also tliose works where the 2;-dependence of the system is averaged out for a 
thin disk[8]. Thus, the possibility of presenting a self-consistent MHD description 
for an accreting disk configuration, calls significant attention, especially in view of 
the implication that the disk morphology can have on the formation of jets from 
compact astrophysical structures [9]. For a discussion which properly traces the 
way to face this relevant configuration problem, see [10]. 

This paper presents a solution to the two-dimensional MHD structure of a 
thin accretion disk, in the linear regime, when the role of a Nernst coefficient (the 
so-called plasma thermal power) is taken into account. In particular, we are recon- 
cihng the oscillating behaviour outlined by B. Coppi, with the accretion features 
of the disk. In fact, in our model, radial and vertical velocities are included in the 
problem, still retaining the poloidal character of the currents living in the plasma. 
The resulting equilibrium configuration is associated with a radial velocity which, 
following the poloidal currents, oscillates as a function of the radial coordinates 
and decays in the vertical direction. Despite its oscillating-like character, this local 
radial velocity is responsible for a non-zero accretion rate and it describes a mech- 
anism for the appearance of a steady matter flux, able to enforce the ring-profile 
associated with the crystalline structure. Nonetheless, the oscillating character of 
this radial velocity prevents a real matter in-fall toward the centre of symmetry. 
Thus, a real accretion of the central object seems to require, in this scheme, the 
presence of non-linear MHD features; but the local character of the model does 
not allow to properly characterise the global accretion profile of the central object, 
that could imply a connection with the boundary layer physics[ll]. 

The structure of the paper is as follows. In Section 2 we present a description 
of the main features characterising an accretion disk model. The tools needed to 
construct the configuration scheme are provided. Section 3 is devoted to fix the 
local nature of the model, properly defining the addressed approximations. In Sec- 
tion 4, we build up the system of partial differential equations which governs the 
two-dimensional MHD equilibrium, by assigning the radial, vertical and azimuthal 
force balance. Section 5 is concerned with the analysis of the linearised equa- 
tions, based on the assumption that the magnetic field induced into the plasma, 
by the presence of the toroidal currents is much less in strength than the external 
field, given by the central object. Section 6 faces the discussion of the electron 
force balance equation, by introducing a non-zero and constant Nernst coefficient 
to preserve the presence of radial and vertical non-vanishing components of the 
matter velocity. Finally wc develop some phenomenological consideration on the 
model, mainly aimed to estimate the strength of the radial velocity and the associ- 
ated accretion rate. Furthermore, we fix the condition to neglect the temperature 
gradients into the radial equilibrium. Brief concluding remarks follow in Section 
7. 
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2 Model for an Accretion Disk 



Let us fix here the basic statements concerning the steady MHD regime for the 
specific case of an accretion disk configuration around a compact astrophysical 
object (a mass over the critical value to have a neutron star within a radius of 
few kilometres), which is also strongly magnetised (a dipole-like field of about 
lO^^Gauss). Thus, we deal with a typical pulsar, accreting material from a binary 
companion via a thin disk structure made up of rotating plasma in equihbrium. 

To adapt the MHD equations to the disk symmetry, we introduce cylindri- 
cal coordinates {r,0,2;}. The gravitational potential of the pulsar takes, in this 
coordinates, the form 

X{r,z)^ — / 2 , 2 ' 
V _l_ 

where G denotes the Newton's constant and M the mass of the spherical central 
object, while the self- gravitation of the disk is regarded as a negligible effect. 
To describe the magnetic field, we take the potential vector in the form 

— * 

A = driler + ipe^ + d^Ue^ , (2) 

where Il{r,z) and ip{r,z^) are arbitrary functions, but only the last one is a 
physical degree of freedom, since the magnetic field reads^ 

B — — dzi>er H — drtpCz . (3) 
r r 

The flux function has to be decomposed as ip = ipo + ipi, where ipo accounts 
for the dipole-hke magnetic field of the pulsar (i.e. tpo = fj,r'^{r'^ + with 
/X = const.), while ipi is a contribution due to the toroidal currents rising in the 
disk configuration. Recalling that the axial symmetry prevents any dependence on 
the azimuthal angle of all the quantities involved in the problem, the continuity 
equation, associated to the mass density e and to the velocity field v, takes the 
explicit form 

^dr{revr) + dz{evz)^0, (4) 

which admits the solution 

^ 1 ^ ^ 1 ^ 

ev — — dzQcr + eu)(r , z)reih H — drOcz , (5) 
r r 

B^-component of the magnetic field can also be included in the theory by a generic term 
Ke^, taken, for instanco, as a function of ip (i.e. K = K{ijj)). For tho sako of simplicity and in 
agreement with Coppi (2005) [1] and Coppi- Rousseau (2006) [7], in our treatment we will fix this 
component equal to zero. 
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= 0(r , being a generic function, but for its odd symmetry in the 2;-coordinate. 
This property has to be required in order to ensure that the accretion rate of the 
disk, as averaged over the vertical direction, be non-vanishing, i.e. 

eVrdz = 47re(r , Zo) = 2nl y^O, (6) 

■zo 

where zq is the half depth of the disk and we recall that Vj. < on the vertical 
average, to ensure a real accretion of the disk. 

The different symmetry of the functions ip and G with respect to the z- 
dependencc, prevents to fix a priori a relation between them. The angular velocity 
u!{r ,z^), describing the differential rotation of the disk, is an even function of z 
and, by virtue of the corotation theorem[12], can be taken as a function of the flux 
surfaces, i.e. u = uj{tjj). 

The momentum transfer through the disk structure is characterised by an az- 
imuthal friction between the radial layers, properly described by a viscosity coeffi- 
cient ?7, responsible for the corresponding non-vanishing stress tensor components. 
The coefficient r] can be decomposed as rj = rjo + rji, where the rii component 
is associated with the additional effects arising from the toroidal currents. The 
necessity to include the electromagnetic back-reaction terms in the viscosity coef- 
ficient, is a consequence of the specific coupling such back-reaction estabhshes. In 
fact, it was demonstrated[l] that the currents induced in the plasma of the disk, 
provide an intrinsic local coupling between the radial and the vertical equilibrium, 
absent in the lowest-order approximation. 



3 Local configuration of the thin disk 

Let us now specialise our model to the case of a thin disk, i.e. Zo/r -C 1 over the 
whole configuration. On average, the dominant contribution to the angular velocity 
of the plasma particles, is then the (equatorial) Keplerian value ujk = (GM/r^)^. 

For a thin disk, the 2;-component of the velocity can be properly identified (see 
Shakura (1973) [3]) with the sound velocity Vs = ^/2KBT/mi {Kb denoting the 
Boltzmann constant, T the plasma temperature, and rrii the ion mass), entering 
the fundamental inequahty 

(7) 

r 

which ensures the possibility to neglect the vertical motion with respect to the 
azimuthal one. 

To estimate on average the ratio between the radial velocity and the azimuthal 
one, we can follow the simple equilibrium condition to fix the asymptotic radial 



5 



velocity of a plasma element falling into the disk 



2 Vr OJ 

VcVr ^ 00 r — OJVs =^ — ~ — , (oj 

Vc being a characteristic frequency of particle collision. A rough estimation pro- 
vides Vc ~ Vgn^^^, where n denotes the number density. Taking Vg ~ 10~^c and 
n ~ 10^*^ cm^^, we get Vc ^ lO^Hz, a value much greater than the typical Keplerian 
frequency. Thus, we conclude that v^i) ~ oukt is, in practice, responsible for the 
main matter flow within the disk, despite its stationary accretion takes place along 
the radial direction. 

In this limit of approximation, the MHD condition on the balance of the Lorcnz 
force has dominant radial and vertical components, providing the electric field in 
the form expected by the corotation theorem, i.e. 

-> V -> (i$ -> OJ 

E^--AB^ -^^^ = -- (driver + d.ijje,) , (9) 

$ denoting the electrostatic potential. However, the axial symmetry prevents a 
dependence of $ on the azimuthal angle and hence, the corresponding 0-component 
of the electric field vanishes identically. This fact requires an additional effect to 
be included into the problem, able to balance the non-vanishing Lorenz force in 
the azimuthal direction, due to the velocity components Vr and v^- We will address 
this crucial question in Section 6, in the limit of a linear theory. 

In the thin disk approximation, the strength of the gravitational field G reads 

as 

G — —ouj^rCr — ouj^zCz . (10) 

We see that a thin disk configuration is justified only in the presence of a sufficiently 
high rotation to deal with a confining vertical gravitational force (for the linear 
scheme, here addressed, the Lorenz force can not affect this statement [7]). 

According to the analysis in [1], we now develop the disk configuration around 
a given value of the radial coordinate ro, limiting our attention to a narrow enough 
interval to express the local angular velocity field in the form 

a; ~ (5a; ~ a;^ + -tt-V'i = ojk + <^o'0i , (H) 

where ipi — ipi{ro , z'^ , r — ro). We approximate the dipole surface function, at ro, 
as V'o(ro) — ^. We drop the rehc ^-dependence in view of the thin nature of the 
disk. 

In what follows, the analysis is performed by addressing the drift ordering ap- 
proximation that fixes the dominant character of the second spatial derivatives of 
-01. In particular, the profile of the disk, so outlined, will include viscous features. 
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as in [3], but making account for the crystalline structure derived in [1]. Indeed, 
reconciling the accretion feature of the disk with the plasma effects, as described in 
a bidimensional MHD approach, we will need a significant deviation from the stan- 
dard mechanism of angular momentum transfer, in the sense that the azimuthal 
equation holds at the level of electromagnetic back-reaction only. 



4 Configuration Equations 

In order to fix the profile of the disk around the configuration at Tq, we have to pro- 
vide the equations governing the radial, the vertical and the tangential equilibrium 
in the presence of the toroidal currents. The radial and the vertical configurations 
are not significantly affected by the viscosity, since its presence mainly concerns the 
differential rotation of the disk[3]. Therefore, these systems stand in the same form 
as in [1], while the implications due to a non-zero viscous stress are summarised 
by the azimuthal equihbrium. 

In order to cast the whole system, we split the energy density and the pressure 
in the form 

e = e(ro, z^) + e(ro, z^, r - Tq) , p = p(ro, z^) + p(ro, r - Tq) , (12) 

where the barred quantities are the contributions existing in absence of the toroidal 
currents, while the terms denoted with a hat are induced by such an electromag- 
netic reaction. According to the framework we outlined, the vertical equilibrium 
is governed by the two relations 

D{z') = = exp"% , 6o(ro) ^ e(ro, 0) , H', = , (13) 

eo(roj rriiLuj^ 

d,p + ulze + + dli^i) = (14) 

47rrQ 

The radial equation underlying the equilibrium of the rotating layers of the disk, 
takes the form 

uj ujk + 5uj <^o(V'o) + TT^^i , (15) 

0-00 AtTTq ^ ' 



The azimuthal equation exactly reads 



evrdr{u!r) + evzdz{ujr) -\- ecuvr = -^dr (jjr^drcu) + dz [r]dz{u!r)] . (17) 
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By using equation (5), the tangential equilibrium above easily stands as 

2 

— dzQdrCu + drOdzCu dzQuj — SrjdrUj + r [drrjdrUj + dzTjOzUj + rj (d^cu + d'^uj)] . 

(18) 

Accounting for the corotation theorem cu — uj{'il)) we can restate the spatial deriva- 
tives of c<j, in terms of the corresponding ones taken on U'. Recalling the expression 
^7 = ^0 + ^1, we can now make the reasonable assumption that the viscosity cor- 
rection Tji be written as 

?7i(ro,'0i) ~ ( ^ ) V'i = ?7oV'i, (19) 

Vi=0 ,r=ro 

we recall that for vanishing ipi, the correction rji vanishes too. In the work [3], 
a valuable proposal for the expression describing the viscosity coefficient 770(^0) is 
provided, i.e. 

2 

^o('^o) = -o.e^)VsoZo , (20) 

where Vsq denotes the sound velocity on the equatorial plane and a is a parameter, 
whose value must be assigned. Putting together all these statements, the tangential 
equation rewrites 

ro u' 

ro [9^V'i^o^ro'0o + droVodril^i + v'o {dri^if + v'o (dzi/^if + Voiro) {dr^i + dl^l^i)] .(21) 

Since we have drofjo ~ 770/^0 ) the first term on the right-hand-side is negligible in 
comparison to the term containing second derivatives. Furthermore, the request 
Vi ^ ^0 imphes that tJq ^ ^0/^1 , ensuring that also the quadratic terms in the 
derivatives of •01 are much smaller than the remaining ones. 

Therefore the form taken by the azimuthal equation for the steady state of the 
disk reads as 

ro [dri^ivMo + Voiro) {d'ri^i + d^iP,)] . (22) 

This equation provides a differential relation between the fiux surface ipi and the 
function © characterising the radial and vertical matter fluxes. 

5 Configuration in the linear approximation 

Let us analyse the compatibility between the radial and the azimuthal equations, 
in the limit when the induced 2;-field Bl is much smaUer than the source one Bzo, 
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I.e. 

^~A;oro^<l, (23) 
where the magnetic surface ipi admits the explicit dependence 

V'l = V'l (^ro, ko{r - ro), . (24) 

Here, we set fep = ^^k/'^-'ao^ with v\q = B'^q/A-k€o being the Alfven velocity in 
the plasma. Furthermore, A denotes a narrow interval for the localisation of the 
2;-dependence. Under these hypotheses, the radial equation rewrites 

9,2V'i + 9,Vi = -^o^(-2>i- (25) 

At lowest orders in z, we take the following expansion D{z'^) ~ l — z'^/H^, being 
of the same order of magnitude of Hq. In this hmit, the radial equation admits the 
solution discussed in [1] which oscillates as a sin function along the radial direction, 
while it exponentially decays in z'^ along the vertical configuration, i. e. ipi reads 
as 

-01 = -0? sin [k{r - ro)] exp{ — ^} , (26) 



where we set k = /coy 1 — and we obtained the identification A = Hn/ko. 
Such a structure is the hnear feature of the crystalline morphology induced in the 
disk by the toroidal currents. 

In the approximation where koVo 2> 1, we can require that the contribution 
due to the currents on the viscosity coefficient be sufficiently small to neglect the 
term in rjQ of (22), i. e. by virtue of the inequality 

^« fc.ro. (27) 
Vo 

Thus, the final form we address for the azimuthal equation is as follows 

- dzQ (dr,^o + ^ + driJi] + drOdz^, = ro77o(ro) + dl^^) . (28) 

We use this equation to complete the configuration scheme of our thin disk equilib- 
rium. Under the same conditions, fixed above for the linear regime, the azimuthal 
equation (22) is restated as follows 

- d^e (dr^t/jo + ^] + drQdz^, = 77oro (^.^V'l + ■ (29) 
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Let us now assume that the term containing the derivative (9^0 be negligible (this 
is natural because it is multiplied by the small r-component of the magnetic field) 
and then divide the equation by dr^^tpQ. Observing that u^drQipo = dr^ujo — — 3a;o/2 
(where ujq = ojk), we rewrite (29) as follows 

Comparing the equation above with (25), we arrive to the relation 

(^o{rQ)vr = -—d^Q = 3r]okoY , (31) 

To 

in which we made use of the definition of the dimensionless function Y = /coV'i/^roV'o- 
Comparing the two equations, we neglected the factor D{z'^) with respect to the 
2;-dependence of V'l, because Hn 3> 6, i.e. koHn 1. 

Such an accreting profile provides a net increasing of the disk mass, according 
to the relation 



eVrdz ~ 6V2'K^^'^rjo{koro)^ -r sin k{r — tq) , 

■oo A* 



(32) 



where we made use of the dipole-like expression for -00 and we extended the inte- 
gration over the vertical coordinate up to infinity. 

It is worth noting that, despite the radial velocity oscillates, the radial matter 
flux is, on average, an in-falling one. In fact, averaging the expression above 
between two nodes around ro, we get 

(M,) = 12V2n'/'rio{korof^ . (33) 

Equation (31) can be easily solved for 6 in agreement with the adopted approxi- 
mation of neglecting the term containing 9^0, as foUows 



© = — 377oA;oro / Ydz . 



Hence we get the vertical velocity in the form 



J Ydz . (34) 



v^^—drQ^-^^ I Ydz. (35) 



^0^0 ^0 

We see that, in the linear approximation, the azimuthal equation ensures the 
existence of a non-zero accretion rate, which is modulated by the oscillating proflle 
of the disk configuration. 
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6 The Electron Force Balance Equation 



Despite we have properly argued that the radial and vertical components of the 
fluid velocity arc significantly smaller than the 0-component of the velocity, nev- 
ertheless, the MHD electron force balance, restated here 

E+-AB^O, (36) 
c 

requires {v A i?)^ = 0, because of = 0, by virtue of the axial symmetry. 
For instance, this feature arises naturally if the corotation theorem is addressed. 
More simply, this request comes out from the independence of (6 characterising the 
electrostatic potential. Thus, on the left-hand side of equation (36) an additional 
term has to appear in order to save the model consistence. In the hmit of the hnear 
theory, the most natural term to be included in the analysis, preserving the MHD 
scenario we addressed, is a weak collisional effect, as phcnomcnologically described 
by a non-zero constant Nernst coefficient. However, to introduce this hypothesis 
we need to give up the idea, pursued above, that the plasma in the disk be perfectly 
isothermal. Thus, we now introduce temperature gradients[13] in the spirit that 
they do not affect significantly the three configuration equations developed in the 
previous sections. Wc discuss only the implications for the electron force balance 
when a Nernst coefficient is present and then we fix the restriction on the model 
parameters to get fuU consistence. 

According to this idea, we split the temperature as T = Tq -|- Ti(r, z'^), with 
To = const, and Ti -C Tq. In the presence of a non-zero Nernst coefficient jV", the 
equation (36) takes the form 

E+-AB= MB A VT . (37) 
c 

Assuming that the gradients of the temperature are negligible in the radial and 
vertical component of this equation, its ^-component provides 

Vr = -AfcdrTi , V, = -Afcd.Ti . (38) 

We now have to observe that the temperature of the plasma in the static MHD 
framework, must obey the equation 

v-VT+^TV-v^O. (39) 

Expressing the gradients of the temperature by the relation (38), the equation 
above takes the form 

hirvr) + d^v^ = {v^^ + vl) . (40) 
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It is easy to see that this equation is automatically satisfied in the linear regime, 
by virtue of the expressions for the radial and vertical velocity (31) and (35) 
respectively. 

The complete consistence of our linear model is then guaranteed by requiring 
that the radial temperature gradient be negligible in the corresponding configura- 
tion equation. Such gradients arise in the radial equation from the spatial variation 
of the pressure p = 2eKBT/mi (from which we getp = 2i^B^oe/"^i+2eo-f^sTi/mj). 
The condition for the negligibility of this term, can be stated, by virtue of equations 
(31) and (38), in the form 



9.0^0 »roJ^, (41) 
V niiJv c 

To get this result, we made account for the fact that obtaining equation (25), we 
divided the original radial equation by the term drQipo/rl- By making use of the 
expression (20) for 770, we can rewrite the condition above in terms of the Alfven 

and equatorial sound velocities as va ^ Svso, with S = ■\J~^^^ (we approximated 
everywhere the sound velocity with its value on the equatorial plane and we have 
taken zq = Hq). Furthermore wc implicitly required that the radial variation of e 
be smaller or, at most, equal to the term with the temperature gradient. 
Finally, the vertical equation, in the linear approximation reads as 

dJ + ulze^O, (42) 

which easily rewrites as an equation for e, i. e. 

-Oze + -^ze -\ — / Ydz — 0. (43) 



rrii 



2 m,A/"c 



Such an equation fixes the following form for the perturbed mass density e = 
e{z^)smk{r — Tq), where the 2;-dependencc comes from the ordinary differential 
equation 

KsTode^u^^_^ 3KsVohk ^o /exp = , (44) 



rrii dz 2 rriiMc J \ 

having defined F° = A;oV^?/i9ro^o- The above equation, using the expression (20) 
for 7^0, can be rewritten in the dimensionless form 



>1, 6= ; Y 
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Figure 1: Behaviour of D for several values of and ^2- The figure outlines an 
increasing behaviour of the function D, and a maximum sometimes arises in the 
range of validity for the linear approximation. 

We conclude this analysis, by stressing that, because of (20), the expression 
(31) provides the following estimation for Vr 

Vr ~ 2aVsokoHoOiY'^) . (47) 

We see that the accretion velocity is few orders of magnitude less than the sound 
velocity on the equatorial plane. In this respect, we observe that the parameter a 
can be taken in the range 10~^ — 1, while the dimensionless term fco-f^o ~ koHn ~ 
^/J3 can be even much greater than unity and the linear approximation requires 

r° < 1. 

Despite the picture we fixed above is well-grounded in the hnear approximation, 
it is well-known that the Nernst coefficient is obtained in a kinetic theory via an 
expansion in the inverse powers of the cyclotronic frequency and it corresponds 
to a first order approximation scheme. Thus, this effect is, in principle, small and 
we can expect that it be the proper explanation to the accretion features of a 
crystalline disk, in the weak field hmit only, i.e. we have to require /5 > 1, to 
ensure the model reliability. 
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7 Concluding Remarks 



Our analysis was devoted to a study reconciling the crystalline structure of the 
disk, as outlined in [1], with the existence of an azimuthal equilibrium configura- 
tion. In fact, the impossibility to neglect the radial component of matter velocity, 
requested by the notion of accretion, leads to include in the equihbrium problem 
the angular momentum transport across the disk and the unavoidable viscous fea- 
tures, associated to the differential rotation. The problem to get a self-consistent 
scheme, in which the electron force balance equation properly holds even along 
the azimuthal direction is then addressed. In fact, the appearance of the radial 
and vertical component of the matter velocity, together with the vanishing nature 
of the azimuthal electric field in an axisymmetric configuration, rise the question 
how to get a consistent equilibrium balancing the Lorenz force along the toroidal 
symmetry. The proposal we investigate here concerns the role that the Nernst 
coefficient can have in this perspective. In agreement with the idea that this is a 
low cyclotronic frequency effect, we treat only the linearised MHD configuration. 
The presence of such a collisional term allows to link the radial and vertical ve- 
locity components to the temperature gradients, as soon as the azimuthal electron 
force balance equation stands. The hnear behaviour is completely self-consistent 
and the adiabatic equation for the plasma pressure is satisfied as a consequence 
of the incompressible nature of the linear MHD, out-coming in this axisymmetric 
scenario. Finally we were able to establish the condition ensuring that the radial 
temperature gradient be neghgible in the corresponding equihbrium, while the lin- 
ear vertical equilibrium is solved exactly. The main issue of our analysis is that 
the radial velocity has an oscillating character, with an amplitude few orders of 
magnitude below the sound velocity. Indeed, we deal with a non-zero accretion 
rate because the material outside ro is greater than the one inside (we are estimat- 
ing the matter density with its value at ro). Such an oscillating profile is a direct 
consequence of the hnear crystalline structure of the model. More than with a real 
matter in-fall, arising expectably in the non-linear regime, we get a matter fiow 
which explains the formation of the ring configuration within the disk. However 
our goal opens interesting perspective on the understanding about the existence of 
a ring-like decomposition in presence of a mechanism of accretion, i.e. significantly 
non zero radial in-fall of the plasma. 
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